Magnetic Scanning Tunneling Microscopy with a Two- Terminal Non-Magnetic Tip: 

Quantitative Results 
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We report numerical simulation result of a recently proposed {P. Bruno, Phys. Rev. Lett 79, 
4593, (1997)} approach to perform magnetic scanning tunneling microscopy with a two terminal 
non-magnetic tip. It is based upon the spin asymmetry effect of the tunneling current between a 
ferromagnetic surface and a two-terminal non-magnetic tip. The spin asymmetry effect is due to the 
spin-orbit scattering in the tip. The effect can be viewed as a Mott scattering of tunneling electrons 
within the tip. To obtain quantitative results we perform numerical simulation within the single 
band tight binding model, using recursive Green function method and Landauer-Biittiker formula 
for conductance. A new model has been developed to take into account the spin-orbit scattering 
off the impurities within the single-band tight-binding model. We show that the spin-asymmetry 
effect is most prominent when the device is in quasi-ballistic regime and the typical value of spin 
asymmetry is about 5%. 
PACS numbers: 61.16.Ch, 73.40. Gk, 75.30.Pd, 75.60.Ch 



I. INTRODUCTION 

Imaging the magnetic structures of surfaces down to 
the atomic level is a major goal of surface magnetism. 
Magnetic scanning tunneling microscopy(MSTM) pro- 
vides a way to image magnetic domains on surface. In 
the conventional approach magnetic sensitivity of tunnel- 
ing current has been based upon the spin- valve effect |Q , 
the tunneling current between two ferro-magnets sepa- 
rated by a tunnel barrier depend on the relative orienta- 
tion of the magnetizations of the ferro-magnets. In this 
approach a magnetic tip has to be used. The experimen- 
tal realization of magnetic scanning tunneling microscopy 
based on spin-valve effect was realized by Wiesendanger 
et. al. Q, who investigated a Cr(001) surface with a 
ferromagnetic Cr02 tip, their observation confirmed the 
model of topological antiferromagnetism between ferro- 
magnetic terraces separated by monoatomic steps. They 
measured a spin asymmetry of the order of 20%. Re- 
cently this method have been used to image magnetic do- 
mains |^-|^]. It was shown that, by periodically changing 
the magnetization of tip, it is possible to separate spin- 
dependent tunnel current from the topographic depen- 
dent current and hence the magnetic structure of surface 
can be recorded. Using this method Wulfhekel et. al. || 
studied magnetic domain structure on single crystalline 
Co(0001) surface and polycrystalline Ni surface. In Refs. 
|0-p[, a two dimensional anti- ferromagnetic structure of 
Mn atoms on tungsten(llO) surface was investigated. It 
was shown that the spin-polarized tunneling current is 
sensitive to the magnetic superstructure, and not to the 
chemical unit cell . 

However the MSTM with a magnetic tip has a draw- 
back that the magnetostatic interaction between the tip 
and magnetic sample can not be avoided, which are likely 
to influence the domain structure. In view of this an 
alternative approach was recently proposed to perform 



the magnetic scanning tunneling microscopy with a two 
terminal non-magnetic tip [j. It is based upon Mott's 
spin-asymmetry effect in scattering caused by disorder 
^| . It was shown that due to spin-orbit coupling the tun- 
nel conductance between the ferromagnetic surface and 
one of the tip terminal depends on the orientation of 
magnetization. Because of the spin-orbit interaction the 
intensity of scattered beam depends on the orientation 
of spin-polarization axis of the incidents electrons, i.e., it 
is sensitive to the spin component perpendicular to the 
scattering plane. In other words tunnel conductance is 
spin asymmetric. However to observe this spin asymme- 
try effect, caused by Mott scattering, a three terminal de- 
vice is a prerequisite. Due to the Casimir-Onsager sym- 
metry relation the conductance of a two-terminal device 
has to be symmetric with respect to magnetic field (in 
our case spin plays the role of magnetic field since as far 
as time reversal properties are concerned 'spin' and 'mag- 
netic field are equivalent), this is a requirement imposed 
by the underlying microscopic time reversal symmetry. 
However in case of three terminal device, there is no such 
restriction on the conductance rather a more generalized 
symmetry relation exists involving all terminals as shown 
by Biittiker Hence to perform magnetization sensi- 
tive scanning tunneling microscopy with a non-magnetic 
tip, it is necessary to use a two terminal tip Q. 

In this work we report numerical simulation result 
of the three terminal STM device within the single- 
band tight-binding model, using recursive Green function 
method and Landauer-Biittiker formula for conductance 
fllPf . We have developed a new model to take into account 
the spin-scattering within the single-band tight-binding 
model. 

The paper is organized as follows, in the next section 
we introduce the single-band tight-binding model includ- 
ing the spin-orbit interaction and the three terminal STM 
device. Section 3 briefly describes the method of calcula- 
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tion and in section 4 we presents some numerical results 
and discussion. 



2. MODEL AND METHOD 

A cross section of the system in xy plane, for the cal- 
culation of spin sensitivity of the proposed two terminal 
non-magnetic tip is shown in Fig.l. The system consist 
of three regions, (i) the ferromagnetic lead (labelled 1 
in Fig.l) , (ii) the central region and, (Hi) the two non- 
magnetic terminals (labelled as 2 and 3 in Fig.l). The 
central region is composed of an insulating tip, such as 
those routinely used to perform atomic force microscopy, 
coated on two opposite faces by a thin metallic film. The 
metallic coating has thickness d. This is shown in cen- 
tral region where the empty circles depict vacuum, black 
circles corresponds to insulating sites and the rest corre- 
sponds to metallic sites(hatched circles) and the impu- 
rities (stars). Between the ferromagnetic surface (gray 
circles in Fig.l) and the tip there is a vacuum layer of 
one lattice spacing (empty circles in Fig.l). The tip is 
placed symmetrical with respect to xz plane. Current 
flows along the two faces of tip which makes an angle of 
±45 degrees with the a;-axis. The structure shown in 
Fig.l consists of three semi-infinite leads (— oo < i < 1 
and N x + 1 < i < oo) separated by tip region 1 < i < N x . 
The thickness of metallic coating on the tip is da where 
a being the lattice constant and the cross section of the 
system is (N y a x N z a), where N y and N z are number 
of sites along y and z axis. For numerical calculation we 
have taken N y =N z — 20 , ^=10 and the metallic coating 
on the tip has a thickness of 4 lattice spacing , i.e.,d=4 
as shown in Fig.l. 

We model the system shown in Fig.l as a single-band 
tight-binding Hamiltonian with nearest neighbor hopping 
parameter t. To obtain appropriate form of single-band 
tight-binding Hamiltonian including spin-orbit interac- 
tion, we discretize the following single-band Hamiltonian 
in continuum, on a simple-cubic lattice, 

H = £ + y(r) + f^ (r) ' S + A(W W X °° ' P > (1) 

where the first two terms are usual kinetic and poten- 
tial energies while the third and forth terms represent 
exchange and spin-orbit interaction, respectively, m* is 
the effective mass of electron, A the exchange splitting, 
\x a unit vector in the direction of magnetization of FMs 
and is given by (cos0sin6>, sin0sin#, cos#), a the Pauli 
operator and P the momentum operator. The discretize 
form of Hamiltonian reads, 



<r,r'>a 



where H so is expressed as 

O-aa'Eijk- (3) 

Here c+ CT the creation operator of an electron with spin 
a at site r, e T on-site energy and Ae r+7akjr+; , aj = 
£r+ 7 a k — £r+iyaj i a i is the lattice basis vector along axis 
i, , c^, , denotes the Pauli matrix elements, a so is di- 
mensionless spin-orbit parameter. The dummy indice 
wy takes the values ±.The summation < r, r' > runs 
over nearest neighbor sites. The symbol eyt is the Levi- 
civita's totally antisymmetric tensor, where ijk label the 
three coordinate axis. 

The tight-binding parameters in equation (2) and (3) 
are related with the parameters in equation (1) in the 
following way, 



(4) 



2m*a 2 ' 



(5) 



(2) 



The above tight binding model includes two factors: 
spin-dependent band structure and spin-independent dis- 
order. The band structure takes into account of the 
difference in the density of states and the Fermi veloc- 
ity between the two spin component in the ferromag- 
net. The disorder represents the structural defects in the 
real STM tip and is source of spin-orbit scattering and 
it takes the form of spin independent random variation 
in the atomic on-site energies. In presence of disorder, 
spin-orbit coupling term causes hopping along the diago- 
nal and is the source of spin-flip scattering. In this sense 
this model is equivalent to next-nearcst-neighbor (nnb) 
tight-binding model, except that in the usual nnb tight- 
binding model, hopping amplitude to the next nearest 
neighbor is fixed while in our model it depends on disor- 
der strength and the spin of electron. Hence within this 
model spin-relaxation length is determined by disorder 
strength. 



3. THEORY 



As shown in Fig.l, the ferromagnet, the left face of 
the tip and the right face of the tip is connected to three 
reservoirs at chemical potentials Vi, V2 and V3 respec- 
tively. Let 7i, I2 and -Z3 be the corresponding incoming 
currents in the three terminals HJlGl. 
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FIG. 1. Cross section of tip-geometry shown in xy plane. 
Gray circles denote Ferromagnetic sample, the empty circles 
depict vacuum, black circles corresponds to insulating sites 
and the rest corresponds to metallic sites(circles with down- 
ward diagonal lines) and the impurities (empty stars) in the 
tip. 



The currents are related to potentials by 

q¥=p 



(6) 



The above expression is gauge invariant and the currents 
conservation law J2i h — requires that G pq =G qp be 
satisfied. 

The calculation of the conductance of the structure 
is based upon the non-equilibrium Green's function for- 
malism flriyi^t . When applied to multi- terminal ballistic 
mesoscopic conductor we obtain the following result for 
the conductance Eof 



G 



.„, — — Tr[T q G TpG ]. 



(7) 



Here p and q enumerates the three terminals and the 
upper indices R and A refer to the retarded and advanced 
Green function of whole structure taking leads into ac- 
count. Here T p / q ) self-energy function for the isolated 
ideal leads and are given by T p ^=t 2 A p ( q ), where A p ^ 
is the spectral density in the respective lead when it is de- 
coupled from the structure. The trace is over space and 
spin degrees of freedom, and all the matrices in equation 
(4) are of size (2N y x N z ,2N y x N z ), where N y and N z 
are number of sites along y and z direction and the factor 
2 takes into account the spin degree of freedom. All the 
quantities in the above equations are evaluated at the 
Fermi energy. To calculate the required Green function 
we use the well known recursive Green function method 



4. RESULTS AND DISCUSSION 

In this section we present numerical results for a sys- 
tem of cross section (20 x 20) in yz plane and a length 
of 10 lattice spacing along the x direction. Number of 
metallic layers on the tip, i.e., d in Fig.l is taken to be 
4 lattice spacings. The hopping parameter, t ,is same for 
all pairs and set to —1 for numerical calculation. The 
on-site energies in the leads and on the metal coating on 
the tip is set to be zero, while in the vacuum layer it 
is £ V ac = 4.0|t|, and in the insulating region in the tip 
it is tins = 10.0 <|. The Fermi level throughout the cal- 
culation is kept fixed at e/ = 3|£| above the bottom of 
the band. For disorder we consider Anderson Model in 
which a random on-site energy, characterized by square 
distribution of width W, is added to the on-site energy 
of perfect case. In our case disorder is added only in the 
metallic coating on the tip; everywhere else the system 
is perfect. 

Before we go over to the discussion of our results we 
briefly mention the correspondence between the physi- 
cal parameters and the model parameters. The relevant 
physical parameters are mean-free-path, spin-relaxation 
length Fermi Energy and the spin polarization of the fer- 
romagnet at the Fermi level, and the model parameters, 
are on-site energy, hopping energy, exchange splitting 
and spin-orbit coupling parameter. Physical parameters 
are related to the model parameters in the following way, 



t rn — 




1 



t\N 3D (e f )((V-Vr) c 



3£ m \t\ 



2|a so |e/ ' 



P = 



Nl(e f )+Nl(e f ) ^J+A + ^f^A' 



(8) 



(9) 



(10) 



where l m , l so and P are elastic mean free path, spin re- 
laxation length and spin polarization of the ferromagnet 
respectively. Here a is lattice spacing and (...) c represents 
the configuration averaging, other symbols have the same 
meaning as defined in the section 2. Below we present 
some numerical results for one typical realization of the 
disorder , we have not performed disorder averaging. 

In Fig. 2 we have plotted the conductance Gn and G13 
as a function of magnetization angle 9 with respect to 
z axis. We rotate the magnetization in yz plane such 
that the magnetization is always perpendicular to the x- 
axis or in other words the angle (f> does not change and 
has fixed value of 90 degrees. To be specific, when 9=0 
and <fi = 90 degrees, magnetization is parallel to z-axis 
while for 6>=90 and <fi=90 degrees the magnetization is 
parallel to y-axis. We have taken e/ = 3|t|, a so = 0.02 
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and A = 2A\t\ and the Anderson disorder strength is 
W = 11*1. 
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FIG. 2. Fig. 2. Conductance (G12 and G13) versus 6 plot 
for the two terminals for a fixed value of <^=90. The other 
parameters chosen for this figure are e/ = 3 . 1 i | , A = 2.4|i|, 
a so =0.02, and W = l|i|. 

This set of parameters corresponds to a mean free path 
of l m = 80a, spin relaxation length of l so — 25l m and 
polarization is P = 50%. We notice that the conduc- 
tance shows approximately cos(#) as a function of angle, 
which is expected since in our geometry the tip is placed 
symmetrically to xz plane. However because of disorder 
the effective axis in the system does not coincide with 
the chosen spin quantization axis, i.e., z axis, and also 
since the structure considered is three dimensional so the 
scattering plane is not fixed hence the conductance vari- 
ation with magnetization angle does not show an exact 
cosine behavior. Also we notice that the variation of G13 
is opposite to that G12. This is in agreement with the 
underlying microscopic timc-rcvcrsiblc symmetry which 
requires that the two terminal conductance should be 
symmetric under time reversal. 
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FIG. 3. Plot of the average two-terminal conductance 
(G12 + Gi3)/2 versus 6 . The other parameters are same 
as those in Fig. 2. 



To verify this point in Fig. 3 we have plotted, the sum of 
(G12 + Giz)/2, we see that the two terminal conductance 
is symmetric with respect to magnetization angle 9 also 
the magnitude of oscillation is much smaller than cither 
G12 or G13. This is due to the fact that the variation 
of G12 or G13 with 9 is of first-order with respect to the 
spin-orbit coupling, whereas the variation of G12 + G13 is 
of second order with respect to spin-orbit coupling. Actu- 
ally the latter can be viewed as related to the anisotropic 
magnetoresistance of ferromagnets, whereas the former 
is related to the extraordinary Hall effect. 

In Fig. 4 we plot spin asymmetry as a function of po- 
larization of ferromagnet for terminal 2. We have defined 
the spin-asymmetry as, 



Gmax r*m:in 
12 — ( - T 12 



(11) 
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r 12 



where to find G™ 2 and G^ 1 ™ we generate a curve as 
shown in Fig. 3 for each set of parameters and from those 
points we get the corresponding maximum and minimum 
values. This is necessary since the variation of conduc- 
tance with magnetization angle does not follow exact co- 
sine behavior, the maxima and minima need not to occur 
exactly at zero and 7r respectively, we have fixed Fermi 
energy at e/ = 3|i| and a so = 0.02. Different curves in 
Fig. 4 corresponds to disorder strengths W = l\t\ (solid 
line), W = 2 \t\ (dotted line), W = 2.5 |i| ( dot-dashed 
line) and W = 4\t\ (dashed line) corresponding mean- 
free-paths are respectively 80a, 10a, 6a, and 3a. Al- 
though all these curves corresponds to different mean- 
free-paths, however the ratio l so /l m is same for all the 
curves and is equal to 25, since this ratio is determined 
by Fermi energy and spin-orbit coupling strength, which 
are kept fixed here. 
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FIG. 4. Spin asymmetry A, as a function of polariza- 
tion for different disorder strength. The other parameters are 
e f = 3.0\t\ and a so =0.02. 

We see that for a fixed disorder strength the spin asym- 
metry increases linearly with the polarization or in other 
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words spin asymmetry is directly proportional to the po- 
larization of ferromagnet. However for a fixed polariza- 
tion value, spin asymmetry shows a non-monotonic be- 
havior. As we increase disorder strength, spin asymme- 
try first increase and then starts decreasing. This shows 
that the spin asymmetry is maximum when the system 
is in quasi ballistic regime, since the multiple scattering 
destroys the spin asymmetry effect. Which is clearly vis- 
ible in the Fig. 4 where spin-asymmetry is maximum , 
for a fixed value of polarization, at a disorder strength 
W = 2\t\ , corresponding to a mean-free-path of 10a, 
while it is minimum for W = 4\t\ corresponding to a 
mean-free-path of 3a lattice spacings. The order of mag- 
nitude of spin- asymmetry is 5%, which is in good agree- 
ment with the prediction in Ref. 0]. 

In Fig. 5 we have studied the behavior of spin asymme- 
try as a function of spin-orbit coupling parameter a so . 
The other parameters are same as in Fig. 4. We no- 
tice that the spin-asymmetry shows a linear behavior for 
small values of a so < 0.03, for larger a so the linear be- 
havior is no longer seen, because for a fixed disorder, 
i.e., fixed l m , as we increase a so , correspondingly , l so , 
i.e. spin-relaxation path decreases, hence the higher or- 
der effect in spin-orbit coupling starts dominating so we 
no longer observe a linear behavior. Also we see that for 
a fixed a so , spin-asymmetry shows a maximum at a dis- 
order strength of around W — 2\t\. This is in harmony 
with the results presented in Fig. 4. The typical value of 
spin asymmetry is of the order of 5%. So from the results 
of Fig. 4 and Fig. 5 we can say with confidence that the ef- 
ficency of the proposed three terminal STM device would 
be maximum when the device operates in quasi-ballistic 
regime. 
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FIG. 5. Spin asymmetry A, as a function of l m /lso for dif- 
ferent disorder strength. The other parameters are e/ = 3.0|t| 
and A = 2.4|t| corresponding value of polarization is 50%. 

In summary we have developed a new model to take 
into account spin- orbit scattering within the single-band 
tight-binding model. Using this model we have done nu- 
merical calculation of magnetic scanning tunneling mi- 
croscopy with a non-magnetic tip. The order of mag- 



nitude of the spin-asymmetry is about 5%, which is in 
good agreement with the qualitative estimate given in 
0, and the effect is maximum when the device oper- 
ates in the quasi-ballistic regime. The spin- asymmetry 
of the present effect is smaller than the one obtained in 
the spin- valve tunneling structures. However, it has some 
advantages. In particular since the tip is non-magnetic, 
it is insensitive to an external magnetic field . This allows 
one to study the domain structure as function of applied 
field. Furthermore, the problem of the magnetostatic 
interaction between the tip and the magnetic sample is 
avoided, which incase of a magnetic tip would give rise 
to undesirable magnetic forces between the tip and the 
sample and are likely to influence the domain structure. 
Another important advantages of this technique is that 
by measuring separately the currents I2 and J3 of the 
two tip terminals, and by combining them appropriately, 
one can separate the weak magnetic contrast from the 
dominant topographic contrast: the sum I2 + 13 depends 
only on the topography, whereas the magnetic informa- 
tion is contained in the difference I2 + I3. Besides all 
these advantages it has an intrinsic limitation that the 
only in-plane components can be studied and also since 
multiple scattering diminishes the spin- asymmetry effect, 
it is necessary that the device operates in a quasi-ballistic 
regime. However to construct such a tip would be exper- 
imentally challenging. 
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